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Purpose: For Monte Carlo simulation of radiotherapy, x-ray CT number of every system needs 
to be calibrated and converted to mass density and elemental composition. This study aims to 
formulate material properties of body tissues for practical two-step conversion from CT number. 

Methods: We used the latest compilation on body tissues that constitute reference adult male and 
female. We formulated the relations among mass, electron, and elemental densities into polylines 
to connect representative tissues, for which we took mass-weighted mean for the tissues in limited 
density regions. We compared the polyline functions of mass density with a bi-line for electron 
density and broken lines for elemental densities, which were derived from preceding studies. 

Results: There was generally high correlation between mass density and the other densities except 
of C, N, and O for light spongiosa tissues occupying 1% of body mass. The polylines fitted to the 
dominant tissues and were generally consistent with the bi-line and the broken lines. 

Conclusions: We have formulated the invariant relations between mass and electron densities and 
from mass to elemental densities for body tissues. The formulation enables Monte Carlo simulation 
in treatment planning practice without additional burden with CT-number calibration. 


I. INTRODUCTION 


II. MATERIALS AND METHODS 


In treatment planning (TP) of radiotherapy, a patient 
body is commonly modeled to be H 2 0 of variable effec¬ 
tive density, which corresponds to electron density for 
high-energy photons or stopping-power ratio for charged 
particles. The effective density is normally converted 
from computed-tomography (CT) number for attenua¬ 
tion of kilovoltage (kV) x ray. The standard approach 
to construct such conversion functions involves experi¬ 
mental modeling of the x ray and stoichiometric analysis 
of body tissues. 1,2 The effective density is, though practi¬ 
cally successful, 3,4 an approximate concept for radiations 
undergoing complex interactions. 

With the advancement of computer technology, Monte 
Carlo (MC) simulation of radiotherapy is becoming fea¬ 
sible, where a radiation is handled as a collection of 
particles individually interacting with matter of known 
composition according to the basic laws of physics. For 
the modeling of body tissues, Schneider, Bortfeld, and 
Schlegel (SBS) applied the stoichiometric calibration to 
construct functions to convert CT number to mass den¬ 
sity and elemental weights. 5 Their conversion functions, 
though commonly used for research, 6,7 are not applicable 
to the other CT systems. Calibration and maintenance 
of the complex one-to-many relations may prevent MC 
simulation from applying to TP practice. 

For patient dose calculation, TP systems commonly 
use a selectable function to convert CT number to ef¬ 
fective density of interest. Recently, a practical two- 
step approach was proposed for TP with proton and ion 
beams, 8,9 where CT number is only converted to electron 
density that is automatically converted to interaction- 
specific effective densities using invariant relations. In 
the present study, we extend the two-step approach to 
promote MC-based TP practice. 


A. Standard tissues and material properties 


We used the standard body tissue data in ICRP Pub¬ 
lication 110, 10 which is the latest compilation of the 
kind. In the publication are mass density p and elemental 
weights w for 53 standard tissues that fully comprise 141 
organs of Reference Male and Female, for which we de¬ 
rived electron density p e and mass fraction m per person 
(50% male/50% female). 

Ignoring tiny mass of air, below 0.90 g/cm 3 was only a 
lung at 0.384 g/cm 3 occupying 1.4% of body mass, which 
we chose for elemental composition of air-containing tis¬ 
sues. In the 0.90-1.00 g/cm 3 region were an adipose tis¬ 
sue at 0.95 g/cm 3 (33.9%) and medullary cavities includ¬ 
ing bone marrow all at 0.98 g/cm 3 (0.4%). In the 1.00- 
1.07 g/cm 3 region were a muscle at 1.05 g/cm 3 (34.5%), 
spongiosa tissues (0.3%), and many other general organs 
(11.7%). In the 1.07-1.101 g/cm 3 region were a skin 
(4.9%), a cartilage (0.6%), and spongiosa tissues (0.7%), 
which are miscellaneously epithelium, connective, and 
fatty bone tissues. In the 1.101 1.25 g/cm 3 region were 
only spongiosa tissues (5.7%). Above 1.25 g/cm 3 were a 
mineral bone at 1.92 g/cm 3 (5.7%) and a tooth at 2.75 
g/cm 3 (0.1%). 

Among the material properties, we adopted mass den¬ 
sity p as the independent variable and examined its cor¬ 
relation to the dependent variables: electron density p e 
and elemental densities pw of six major elements M = 
{H, C, N, O, P, Ca}. To define the regional representative 
tissues, we took mass-weighted mean for the set of den¬ 
sities p = ( p,p e ,pwH,pwc,pw N ,pwchPWp,pwca ,) over 
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tissue t in region R by 
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For the standard tissues and the regional representative 
tissues, we calculated the residual weight and the mean 
residual atomic number by adding weights and averaging 
atomic number Z r over residual element r, 


Wnes = ^2 Wr and = r ( 2 ) 

with which the residual mass could be approximately in¬ 
cluded in MC simulation. 


B. Density segmentation and tissue mixing 


Compared to muscle/organ tissues, adipose/marrow 
tissues have high concentration of fat. Teeth have high 
concentration of minerals, connective tissues have high 
concentration of collagen, and bones are in between 
them. The concentration varies among and within indi¬ 
vidual tissues and are generally correlated with density. 
Anatomically, adipose tissues are neighboring to muscles 
and organs, muscles are connected to bones via connec¬ 
tive tissues, and teeth are connected to jaw bones. At 
these interfaces, tissue mixing may occur in an image of 
finite spatial resolution. Therefore, we modeled an arbi¬ 
trary tissue of mass density p to be mixture of represen¬ 
tative tissues 1 and 2 that comprise a density segment of 
pi < p < p 2 - The other densities are interpolated from 
those of the representative tissues by 


P = 


P2 ~ P 
P2 - Pi 


Pi + 


P~ Pi 
P2 ~ Pi 
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( 3 ) 


in a mass-weighting manner. 5,11 In other words, we as¬ 
signed the representative tissues to the polyline points 
for conversion from mass density to the other densities. 

For fatty tissues lighter than the representative adi¬ 
pose/marrow tissue, we extended the soft-tissue region 
down to 0.90 g/cm 3 , which is the mass density of human 
fat at 37°C. 12 Similarly, we extended the region for air- 
containing tissues up to 0.80 g/cm 3 , leaving a transition 
segment of 0.80-0.90 g/cm 3 to avoid discontinuity of the 
material properties and to cope with body fat that escape 
over the segment boundary as viewed in CT image. We 
also extended the hard-tissue region up to 3.00 g/cm 3 
for teeth heavier than their representative. To define 
the extended boundaries, we applied Eq. 3 for extrap¬ 
olation although, for the fat, the N, P, and Ca weights 
were forced to zero and the C weight was adjusted to sum 
up to 100% to comply with general composition of fatty 
acid. For the polyline tissues, we compiled mass frac¬ 
tion, mass and electron densities, elemental and residual 
weights, and mean residual atomic number. 



FIG. 1. Correlation between mass and electron densities for 
male (+) and female (O) tissues plotted with the polyline 
function (solid lines), the bi-line function (dotted lines) and 
“adipose 3” (*) at 0.93 g/cm 3 by Hiinemohr et al. (2014) 15 , 
and embedded subplots for the box-shaped areas. 


C. Comparison with preceding formulations 

Hiinemohr et al. formulated a bi-line relation between 
mass and electron densities for the body tissues compiled 
in Refs. 13 and 14. 15 From the same dataset, SBS se¬ 
lected their representative tissues for mixing of air/lung, 
adipose tissue/adrenal gland, small intestine/connective 
tissue, and marrow/cortical bone, 5 from which we de¬ 
rived their relations between mass density and elemental 
densities as broken-line functions. We compared our re¬ 
sults with those preceding formulations. 


III. RESULTS 

Figure 1 shows the correlation between mass and elec¬ 
tron densities, where p ew = 3.343 x 10 23 /cm 3 is the elec¬ 
tron density of water. The two fitting functions, the poly¬ 
line and the bi-line by Hiinemohr et al., were generally 
consistent with the standard tissues except around the fat 
at 0.90 g/cm 3 and around the tooth at 2.75 g/cm 3 , for 
which the bi-line function gave 0.892 (—1.7%) and 2.486 
(—1.3%), respectively. The discontinuity in the bi-line at 
their “adipose 3” was settled by the polyline with bet¬ 
ter fitting. Figure 2 shows the correlation between mass 
density and elemental densities. The polyline functions 
and the broken-line functions according to SBS were both 
generally consistent with the standard tissues except for 
the C, N, and O densities in the 1.0-1.1 g/cm 3 region. 
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TABLE I. Material properties of the polyline tissues: mass fraction m, mass density p, relative electron density p e //Oe w , 
elemental and residual weights W{H,c,N,o,p,Ca} and w Tes , and mean residual atomic number Z Tea , based on ISO Standard 2533 
(1975) 16 for the air and on ICRP Publication 110 (2009) 10 for the regional representative tissues. 


X® Tissue type 

m 

% 

P 

Pe 

WH 

wc 

WN 

wo 

w p 

WCa 

Wres 

Zres 

g/cm 3 

Pew 

% 

% 

% 

% 

% 

% 

% 

1 Air 


0.001 

0.001 

0.00 

0.01 

75.52 

23.17 

0.00 

0.00 

1.30 

18.0 

2 Lung 

1.4 

0.384 

0.380 

10.3 

10.7 

3.2 

74.6 

0.2 

0.0 

1.0 

15.9 

3 Air-containing 


0.80 

0.793 

10.3 

10.7 

3.2 

74.6 

0.2 

0.0 

1.0 

15.9 

4 Fat 


0.90 

0.907 

12.09 

84.39 

0.00 

3.52 

0.00 

0.00 

0.00 

n/a 

5 Adipose/marrow 

34.3 

0.950 

0.952 

11.40 

58.92 

0.74 

28.64 

0.00 

0.00 

0.30 

14.7 

6 Muscle/organ 

46.5 

1.049 

1.040 

10.25 

14.58 

3.20 

70.87 

0.21 

0.02 

0.87 

16.8 

7 Miscellaneous 

6.3 

1.090 

1.077 

9.94 

20.90 

3.84 

63.73 

0.45 

0.27 

0.87 

15.5 

8 Spongiosa 

5.7 

1.137 

1.116 

9.30 

39.15 

2.22 

41.71 

2.36 

4.60 

0.66 

14.9 

9 Mineral bone 

5.7 

1.92 

1.784 

3.6 

15.9 

4.2 

44.8 

9.4 

21.3 

0.8 

13.1 

10 Tooth 

0.1 

2.75 

2.518 

2.2 

9.5 

2.9 

42.1 

13.7 

28.9 

0.7 

12.0 

11 Extra tooth 


3.00 

2.739 

1.93 

8.27 

2.65 

41.58 

14.53 

30.37 

0.67 

11.8 


As shown in Fig. 3, the highest and lowest mean residual 
atomic numbers were 20.6 at 1.04 g/crn 3 for the thyroid 
and 12.0 at 2.75 g/crn 3 for the tooth. The global mean of 
15.95 approximately corresponded to element S. Table I 
shows the resultant material properties of the polyline 
tissues. 


IV. DISCUSSION 

A concise yet complete dataset of the polyline tissues 
may be suited to stoichiometric analysis for CT-number 
calibration. 1,2 In the cases where electron density is avail¬ 
able, the most likely mass and elemental densities can be 
determined by Eq. 3 with variables (p, pi , p%) replaced to 
(p e , p el , p e2 )- The resultant mass density and composi¬ 
tion including residual weight for element S will consti¬ 
tute a volumetric patient model for MC simulation. 

The general agreement between the formulations indi¬ 
cates the consistency between the compilations of body- 
tissue data. The mass-weighting approach of this study 
took advantage of the publication that focused on the 
computational phantoms of reference adult male and fe¬ 
male rather than on variations in age, physical status, or 


individual. 10 

The poor fitting in C, N, and O densities was due to 
undifferentiated inclusion of spongiosa tissues occupying 
1% of body mass in the 10. 1.1 g/crn 3 region. They could 
potentially be resolved by anatomical identification or in¬ 
dependent quantitative imaging, in which case either an 
extended spongiosa/mineral-bone segment or a separate 
marrow/spongiosa segment should be applied to them. 

Recently, megavoltage (MV) CT and dual-energy (DE) 
CT have been investigated for direct electron-density 
imaging. 17,18 Hiinemohr et al. formulated the compo¬ 
sition of body tissues as a two-variable function of elec¬ 
tron density and effective atomic number from DECT 
with improved accuracy. 15 The DECT with kV and MV 
x rays would be ideal. 19 Nevertheless, use of MVCT or 
DECT will potentially mitigate metal artifact and beam 
hardening, which limits the accuracy of kVCT. 20,21 

In conclusion, we have formulated the invariant poly¬ 
line relations between mass and electron densities and 
from mass to elemental densities for body tissues. The 
formulation enables MC simulation in TP practice with¬ 
out additional burden with CT-number calibration. 

The author thanks M. Nakao of HIPRAC, Hiroshima 
for private communication on his inspiring research. 
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Mean residual atomic number, Z res 


FIG. 3. (a) Scatter plot of mean residual atomic number 

versus mass density for male (+) and female (O) tissues with 
thickness varied with residual mass and (b) mass-weighted 
histogram of mean residual atomic number. 





